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SUMMARY 

Research  aimed  at  developing  reliable  methods  to  predict  the 
growth  of  a  crack  in  aircraft  built-up  structures  under  realistic 
load  conditions  is  current  being  carried  out  at  the  Institute 
of  Aeronautics  of  the  University  of  Pisa  under  the  three-year 
research  contract  DA-ERO-7 8-G- 107 . 

This  paper  presents  the  results  obtained  in  the  course  of  the 
first  year.  Such  results  concern  two  main  topics,  namely 

-  prediction  of  crack  growth  in  riveted  stiffened  panels  under 
constant  amplitude  loading 

-  prediction  of  crack  growth  in  plain  sheets  under  variable  am¬ 
plitude  loading. 

With  regard  to  the  first  topic,  the  results  concern  mainly 
the  evaluation  of  rivet  flexibility  and  friction  forces  between 
the  stringer  and  sheet  cover.  Such  quantities  have  been  obtained 
by  analysis  of  crack  growth  test  data  by  means  of  a  SKESA  computer 
program  developed  ad  hoc,  which  establishes  them  by  searching 
for  the  best  regression  by  correlating  the  crack  rate  with  the 
stress  intensity  factor. 

As  far  as  the  second  topic  is  concerned,  a  rationale  has  been 
developed  and  implemented  in  a  CADAV  computer  program  to  eval¬ 
uate  existing  methods  for  computing  crack  growth  under  variable 
amplitude  loading.  Attention  has  been  focused  on  four  methods, 
namely:  non-interactive ,  Wheeler,  Bell-Eidinof f  and  Willenborg 
methods. 

Tests  have  been  performed  both  at  constant  and  variable  am¬ 
plitude  loading .  Constant  amplitude  test  data  has  been  used  to 
obtain  reference  K-rate  relationships .  Variable  amplitude  test 
data  obtained  utilizing  a  FALSTAFF  spectrum,  has  been  compared 
with  prediction  by  means  of  the  CADAV  which  classifies  the  meth 
ods  on  the  basis  of  a  statistical  criterion. 
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1  -  INTRODUCTION 


This  paper  presents  the  first  year’s  results  of  an  investigation 
carried  out  at  the  Institute  of  Aeronautics,  Pisa  University, 
under  contract  DA-ER0-78-G-107 . 

The  objective  of  this  investigation  is  the  development  of  a 
reliable  procedure  for  evalua- ing  the  growth  of  a  crack  in  a 
built-up  structure  under  realistic  loading  conditions. 

To  this  end  tests  have  been  conducted  on  stiffened  and 
unstiffened  panels  to  monitor  the  growth  of  a  crack  both  under 
constant  and  variable  amplitude  loads. 

The  test  data  has  been  evaluated  on  the  basis  of  a  deep  theo 
retical  analysis  which  has  allowed  us  to  reach  significant 
results  on  the  following  main  topics: 

-  influence  of  the  fastener  flexibility  and  stringer-sheet  friction 
forces  on  fatigue  crack  propagation  in  stiffened  panels 

-  reliability  of  existing  methods  for  predicting  the  growth  of  a 
crack  under  loading  conditions  in  plain  sheets. 

2  -  BRIEF  STATEMENT  OF  THE  BASIS  OF  THE  RESEARCH. 

The  need  for  investigation  into  crack  growth  in  riveted  built- 
up  structures  under  variable  amplitude  loading,  springs  from  two 
main  causes,  namely  the  uncertainty  surrounding  the  computational 
methods  of  the  stress  intensity  factor  and  the  complexity  of  the 
crack  growth  interaction  effects  caused  by  load  amplitude  varia¬ 
tions  . 

The  uncertainty  surrounding  the  K  values  has  to  be  ascribed 
mainly  to  the  difficulty  of  assessing  the  actual  values  of  the 
mutual  forces  existing  in  the  junction  between  the  stringers  and 
the  sheet  cover,  (junction  forces) . 

For  a  given  structure  geometry,  the  junction  forces  depend  on 
such  quantities  as  the  rivet-hole  deformation  and  the  friction 


2 


between  stringer-skin  contacting  surfaces.  Only  rough  approxima 
tions  are  currently  available  for  such  quantities  and  therefore 
the  stress  intensity  values  computed  are  affected  by  errors  that 
may  be  sufficiently  large  as  far  as  the  crack  length,  the  struc¬ 
ture  geometry  and  load  level  are  concerned. 

As  a  consequence  of  the  above,  the  basic  problem  of  computing 
the  crack  growth  under  constant  amplitude  loading  also  involves 
uncertainty  in  the  case  of  a  stiffened  structures. 

To  overcome  the  problem  knowledge  on  the  junction  forces  must 
be  improved  to  such  an  extent  as  to  compute  the  stress  intensity 
factor  with  an  accuracy  comparable  to  that  reached  in  the  evalua 
tion  of  K  in  plane  sheets. 

As  far  as  the  second  cause  is  concerned,  the  role  of  the  inter 
action  effects  in  fatigue  crack  growth  problems  is  well  documen¬ 
ted,  ref . | 1 | . 

One  such  problem  which  can  be  taken  as  an  example  is  the  im¬ 
plementation  of  reliable  crack  growth  prediction  methods,  one  of 
the  biggest  problems  in  damage  tolerant  structure  design. 

The  research  effort  in  this  area  has  given  rise  to  very  prom¬ 
ising  results  ^nd  certain  important  prediction  methods  are  cur¬ 
rently  available, | 2,3,4, 5 ] . 

However,  the  overall  picture  of  damage  accumulation  during 
crack  growth  under  variable  amplitude  loading  is  still  insuffi¬ 
ciently  clear  and  further  investigation  is  required  to  improve 
the  reliability  of  crack  growth  prediction  methods. 

The  production  of  experimental  data  on  crack  growth  under  load 
sequences  representative  of  airplane  operations  must  be  consider¬ 
ed  of  particular  significance  in  this  respect. 

Further,  as  all  the  existing  methods  derive  their  experimental 
evidence  mainly  from  tests  performed  on  sheet  specimens,  there  is 
a  need  for  crack  growth  data  drawn  from  tests  on  representative 
airplane  structures  under  realistic  loading  conditions. 

In  conclusion,  the  success  of  the  search  for  a  reliable  meth¬ 
od  to  predict  crack  growth  in  stiffened  structures  under  variable 
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amplitude  loading  depends  upon  the  solution  of  two  related  prob 
lems,  namely, 

-  further  insight  into  the  problem  of  the  evaluation  of  junction 
forces  to  produce  a  reliable  K  computational  approach  and 

-  critical  evaluation  of  the  state-of-the-art  crack  growth  pre¬ 
diction  methods  for  variable  amplitude  loading  through  the 
comparison  of  a  set  of  experimental  crack  growth  data  repre¬ 
sentative  of  operational  activity  as  far  as  load  sequences  and 
specimen  configurations  are  concerned. 

In  connection  with  the  first  problem  mentioned  above  significant 
results  can  be  obtained  by  suitable  analysis  of  constant  amplitude 
crack  growth  data. 

The  second  problem  can  be  conveniently  dealt  with  in  two 
successive  stages. 

In  the  first  stage,  by  keeping  the  specimens  to  be  tested  as 
simple  as  possible  (viz.  plane  sheets  for  which  SIF  solutions 
are  well  known) ,  the  reliability  of  existing  prediction  methods 
can  be  assessed  against  crack  growth  data  obtained  from  different 
load  spectra. 

In  a  second  stage,  the  capability  of  the  most  promising  pre¬ 
diction  methods  is  further  exploited  in  relation  to  crack  growth 
data  drawn  from  complex  structures  (built-up  structures)  under 
spectrum  loading  representative  of  the  operational  activity. 

In  the  present  report  the  results  obtained  in  the  first  year 
of  the  research  will  be  presented. 

Such  results  refer  to  the  following  main  topics: 

a1  -  constant  amplitude  fatigue  tests  on  stiffened  panels, 

a2  -  development  and  implementation  of  a  rationale  to  obtain 
information  on  the  junction  forces  (fastener  flexibility 
and  friction  forces)  from  the  crack  growth  data  in  stiffened 
panels, 

b1  -  crack  growth  tests  under  variable  amplitude  loading  in  plain 
sheets, 
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b2  -  development  and  implementation  of  a  rationale  to  evaluate 
the  reliability  of  existing  prediction  methods  of  crack 
growth  under  variable  amplitude  loading. 

3  -  CRACK  GROWTH  IN  RIVETED  STIFFENED  PANELS  UNDER  CONSTANT 
AMPLITUDE  LOADING. 

In  the  present  section,  attention  will  be  focused  on  the 
results  obtained  up  to  now  concerning  the  improvement  of  fatigue 
crack  growth  prediction  techniques  in  riveted  stiffened  struc¬ 
tures  under  constant  amplitude  loading. 

Research  in  this  area  has  been  mainly  devoted  to  assessing 
the  influences  of  fastener  flexibility  and  stiffener  sheet  cover 
friction  forces  on  the  SIF  and  to  develop  and  substantiate  a  ra¬ 
tionale  to  obtain  a  reliable  evaluation  of  such  quantities. 

3,1  -  THE  INFLUENCE  OF  FASTENER  FLEXIBILITY  AND  FRICTION  FORCES 
ON  THE  FATIGUE  CRACK  GROWTH  PHENOMENON. 

As  explained  above,  to  effectively  predict  fatigue  crack  growth 
in  a  built-up  structure,  it  is  necessary  to  evaluate  fastener 
flexibility  and  the  junction  friction  forces. 

A  lot  of  information  on  both  the  above  mentioned  points  can 
be  deduced  by  comparing  the  expected  K-rate  relationship  and  the 
relationship  based  on  crack  growth  test  data  and  K  values  computed 
on  the  basis  of  selected  fastener  behaviour  (elastic, or  elastoplastic 
behaviour  with  a  selected  value  of  flexibility)  and  given  friction 
forces . 

Fig.1  gives  a  qualitative  indication  on  the  way  the  information 
can  be  obtained  from  the  said  comparison.  In  the  upper  part  of  the 
figure  typical  data  obtained  from  a  plain  sheet  test  is  shown. Here 
both  the  crack  rate  da/dn  and  the  stress  intensity  factor  K  are 
increasing  functions  of  the  crack  length  a.  Therefore,  the  rela¬ 
tionship  between  da/dn  and  K  is  biunivocal  and  each  point  of  the 
curve  corresponds  to  an  single  value  of  the  crack  length. 
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In  the  case  of  stiffened  panels,  both  da/dn  and  K  show  maxima 
and  minima  according  to  variations  in  the  crack  length,  so  that 
the  same  point  of  the  da/dn  versus  K  curve  may  correspond  to  more 
than  one  value  of  the  crack  length. 

In  these  conditions,  anomalies  of  the  da/dn-K  curve  are  pos¬ 
sible  due  to  errors  or  in  the  computational  methods  for  K  (cen¬ 
tral  part  of  the  figure)  or  in  the  evaluation  of  da/dn  (lower 
part  of  the  figure) . 

In  all  cases,  the  da/dn-K  curve  is  no  longer  a  single  value 
relationship  and  takes  on  the  different  shapes  shown  in  a  qual¬ 
itative  way  in  the  central  and  lower  part  of  Fig.1. 

These  anomalies  in  the  da/dn-K  curves  can  be  profitably  used 
to  gather  information  on  junction  forces  (fastener  and  friction 
forces) ,  if  a  highly  accurate  procedure  is  available  to  compute 
da/dn.  In  such  conditions  all  the  anomalies  in  da/dn-K  curves 
obtained  by  processing  crack  growth  test  data  can  be  traced  back 
to  errors  in  the  SIF  computational  approach.  Since  the  only 
significant  causes  of  error  in  the  computation  of  the  SIF  can  be 
ascribed  to  an  incorrect  value  selection  of  the  fastener  flexibil 
ity  and  friction  forces,  it  is  possible  in  principle  to  establish 
that  the  "best  values"  of  both  these  quantities  are  those  which 
eliminate  or  minimize  such  anomalies. 

A  first  application  of  this  concept  was  shown  in  references 
6  and  7,  utilizing  a  very  rough  approximation  for  the  computa¬ 
tion  of  the  stress  intensity  factor,  namely  a  zero-friction, 
rigid-fastener  idealization. 

The  anomalies  of  the  da/dn-AK  curves  so  obtained  allowed  us 
to  establish  a  qualitative  classification  of  the  influences  of 
different  parameters  (rivet  types,  stiffener  types  and  materials) 
on  the  fatigue  crack  growth  in  riveted  stiffened  structures. 

3,2  -  FASTENER  FLEXIBILITY  AND  FRICTION  FORCES  EVALUATION  PRO¬ 
CEDURE. 

To  obtain  a  quantitative  assessment  of  these  influences  the 
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following  rationale  (Fig. 2)  has  been  developed.  Crack  growth 
data  is  obtained  by  testing  stiffened  panels  under  constant 
amplitude.  The  da/dn  fatigue  crack  propagation  rate  is  then 
calculated  with  a  spline  function  interpolation  method  (Ap¬ 
pendix  1).  The  stress  intensity  factor  is  computed  using  dis 
placement  compatibility  methods  obtained  by  modifiyng  the 
original  Poe  approach  by  taking  fastener  elastic  flexibility 
and  plasticity  into  account  and  by  allowing  for  the  influence 
of  friction  between  sheet  and  stringer.  Rivet  flexibility  in 
the  elastic  range  is  taken  into  account  by  means  of  a  parameter 
£  which  is  the  ratio  between  the  actual  flexibility  and  a  computed 
value. 

In  the  present  research  the  computed  values  have  been  obtained 
on  the  basis  of  the  simplified  theory  as  summarized  in  Fig. 3  and 
of  the  Swift  formula,  ref.|8|. 

The  fastener  load  which  starts  significant  plastic  deformation 
is  calculated  on  the  basis  of  simple  lap  joint  tests. 

The  friction  forces  have  been  simulated  by  a  set  of  fictitious 
rigid-perfectly  plastic  fasteners  placed  between  the  actual  fast 
eners. 

The  first  fictitious  fastener  is  placed  half  way  between  the 
first  and  the  second  rivet,  the  second  fictitious  fastener  between 
the  second  and  the  third  rivet  and  so  on  as  shown  in  Fig. 4.  The 
fastener  load  Pf,  in  correspondence  with  which  sliding  takes 
place,  is  considered  a  variable  to  be  determined. 

For  each  couple  of  values  £  and  Pf  the  stress  intensity  factor 
is  computed  for  all  the  crack  lengths  relevant  to  the  test  under 
examination. 

The  da/dn  crack  rate  computed  by  means  of  the  spline  approach 
is  correlated  with  AK  obtained  following  the  above  mentioned 
procedure.  A  best-fit  analysis  is  then  carried  out  utilizing  a 
Paris  K-rate  relationship.  The  best-fit  is  then  evaluated  through 
such  statistical  quantities  as  the  standard  deviation,  correlation 
coefficient  and  F-data.  The  best  values  of  £  and  Pf  are  then 
determined  when  extremum  conditions  (minimum  for  standard  devia¬ 
tion  and  maximum  for  the  other  two  quantities)  are  found. Obviously, 


only  one  statistical  quantity  could  be  used  to  obtain  the  best 
values  of  £  and  Pf  but  redundant  information  is  felt  to  be 
useful . 

This  rationale  has  been  implemented  by  means  of  the  SKESA 
computer  program  based  on  three  main  modules,  namely  the  spline 
module,  the  stress  intensity  module,  and  the  best-fit  statistical 
evaluation  module. 

By  means  of  such  a  computer  . t  /gram  the  a-n  test  data  is 
processed  until  the  best  combin_: ion  of  £  and  Pf  is  found. 

3,3  -  experimental  crack  growth  data. 


The  rationale  described  in  the  previous  section  was  applied  by 
means  of  the  SKESA  computer  program  to  a  set  of  a-n  data  obtained 
by  constant  amplitude  crack  propagation  tests. 

The  main  characteristics  of  the  stiffened  panels  are  shown  in 
Tab. I  and  Tab. II.  All  the  panels  are  stiffened  by  straps  of  the 
same  material  as  the  sheet  cover  riveted  with  countersunk  or  round 
head  rivets.  Both  2024-T3  and  7075-T6  aluminum  alloys  were  used  in 
the  panel  construction. 

The  stiffened  panels  of  Tab. I  differ  regarding  the  geometric 
configuration,  the  type  of  rivet  and  the  material.  The  panels  of 
Tab. II  differ  only  as  far  as  the  type  of  rivet  is  concerned. 

The  panels  of  Tab. I  were  tested  in  the  course  of  previous 
research  and  the  relevant  test  data  has  already  been  utilized  for 
preliminary  analys  s  of  the  crack  growth  phenomenon  as  outlined 
in  ref. | 6,7 | . 

The  panels  of  Tab. II  were  tested  in  the  course  of  the  present 
investigation  to  obtain  further  data.  The  tests  were  performed 
with  an  improved  testing  technique  (see  Appendix  2)  based  on  a 
servo  controlled  loading  machine  and  on  measurement  of  the  suc¬ 
cessive  crack  tip  positions  obtained  by  means  of  a  sliding  micro 
scope  which  allowed  us  an  absolute  error  on  the  crack  length  no 
higher  than  .1  mm.  and  pratically  no  error  on  the  number  of  cycles. 
The  results  of  such  tests,  which  are  given  in  Tab. Ill  and  plotted 
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in  figs  5a  and  5b,  can  also  be  profitably  used  to  check  how 
the  accuracy  of  the  crack  length  and  load  measurement  can  in¬ 
fluence  the  scatter  of  the  phenomenon. 

3,4  -  CRACK  GROWTH  DATA  ANALYSIS. 

All  the  a-n  data  relevant  to  the  panels  in  Tab. I  and  II  have 
been  analyzed  utilizing  the  SKESA  computer  program. 

The  elastic  flexibility  of  the  rivets  has  been  evaluated  by 
means  of  a  parameter  given  by: 

dP/dA  *  Erd 

Ci  =  - 

Fi(Xi,ui) 

where  the  first  factor  of  the  numerator  is  the  actual  flexibility 
(see  figs  2,3)  and  the  denominator,  part  of  the  1/Erd  factor, is  a 
reference  value  of  the  flexibility  function  of  the  non  dimensional 
geometric  ratios  A j_  and  the  Young  moduli  ratios  p-^. 

Two  different  types  of  function  have  been  selected  for  the 
purposes  of  the  present  analysis,  namely  the  function  Flf  defined 
in  the  lower  box  of  Fig. 3  which  corresponds  to  ^lfand  the  function 

Er  d  Er  d 

F2=  5  +  0.8  —  —  +  0.8  —  — 

Es  ts  Ec  tc 

which  corresponds  to  £2  obtained  from  ref.  8  and  referred  to  hereafter 
as  the  Douglas  flexibility  formula. 

The  plastic  range  of  the  load-displacement  relationship  is 
accounted  for  by  a  very  slow  slope  linear  variation  starting  from 
a  load  P  which  depends  on  the  shape  and  dimensions  of  the  rivets. 

The  fictitious  rivets  which  simulate  the  friction  forces  have 
been  represented  by  the  sliding  force  Pf  (see  Fig. 2)  by  means  of 
the  parameter 

Pf 

n  =  - 

Sy  Ar 
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where  Sy  is  the  yielding  stress  of  the  rivet  material  and  Ar  the 
rivet  shank  cross  section  area. 

Therefore, n  embodies  the  scatter  due  to  the  variability  both 
in  the  normal  force  and  in  the  friction  coefficient. 

To  illustrate  the  type  of  data  processing  performed  by  means 
of  the  SKESA  computer  program,  some  typical  results  are  shown  in 
figs  6a  to  8b. 

Fig. 6a  shows  the  influence  of  the  friction  forces  on  the  com¬ 
puted  values  of  the  SIF.  Fig. 6b  shows  the  influence  of  flexi¬ 
bility  on  the  same  quantity. 

Fig. 7  describes  how  the  shape  of  the  da/dn-  AK  curve  modifies 
with  varying  £  and  n • 

Each  plot  shows  da/dn  (obtained  by  test  data  and  the  spline 
module)  as  a  function  of  AK  (obtained  by  the  SIF  module  for 
selected  values  of  £  and  n),  together  with  the  best  fit  curve 
computed  forcing  the  da/dn-AK  to  conform  to  a  Paris ^  type  K- 
rate  relationship. 

Lastly,  figs  8a  and  8b  show  typical  plots  of  F-data  and  stand¬ 
ard  deviation, o,  as  a  function  of  £  and  n;  such  plots  are  used  to 
obtain  the  best  values  of  £  and  namely  the  values  £  and  n, which 
maximize  (minimize)  F-data  (the  standard  deviation) . 

To  contain  the  computational  efforts  within  acceptable  limits, 
only  a  relatively  small  number  of  values  of  the  variables  £  and  n 
were  used  to  draw  such  plots. 

The  couple  of  values  £,n  which  gives  extremum  conditions  among 
the  values  selected  for  £  and  n»  was  assumed  to  be  the  best  values 
of  £  and  n  without  trying  a  better  allocation  of  the  extrema  by 
means  of  some  extrapolation-interpolation  technique. 

In  this  way  all  the  values  which  lie  in  a  given  range  around 
the  selected  values  of  £  and  n  are  considered  equivalent. 


(1)  The  Paris  type  K-rate  relationship  has  been  selected  instead 
of  other  types  of  relationships  (Forman  or  Collipriest ) ,  since 
it  works  accurately  in  the  hK  range  typical  of  the  present 
test  data. 
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This  approach  is  felt  to  be  adequate  for  present  purposes 
since  £  and  n  ,  being  random  variables,  can  be  characterized  by 
their  statistic  distributions.  And,  if  one  tries  to  determine 
such  distributions  simply  through  histograms,  he  only  needs  to 
know  how  many  values  of  the  random  variable  lie  within  each 
interval  selected  to  construct  the  histogram. 

On  the  basis  of  such  considerations,  as  far  as  n  is  concerned, 
after  preliminary  trials  the  following  values  were  considered  sufficien 

nj  =  0,  n  2 =  0.07,  n  3=  0.014,  ru,=  0.021 
which  correspond  to  sliding  forces  of 

Pfl  =  0,  Pf2  =  20,  Pf3  =  40,  and  Pf4  =  60  Kg, 

on  the  basis  of  a  rivet  diameter  of  1/8  and  Sy  =  36  Kg/mm2. 

In  this  way  each  value  m  can  be  considered  to  represent  all 
the  values  lying  approximately  in  the  interval 

ni  -  0.035  £  n  +  0.035 

As  far  as  the  selection  of  the  variables  £  is  concerned,  in  a 
first  trial  we  worked  with  three  values  of  £ ,  namely 

?!  =  0.3  £  =  0.7  £  3  =  1 

By  inspection  of  the  diagram  obtained  with  such  values  a  pre¬ 
liminary  rough  allocation  of  the  maximum  of  F-data  was  obtained. 

From  such  an  inspection  it  was  found  that  £  is  scattered  within 
an  interval  ranging  approximately  from  0.1  to  1. 

As  a  consequence  of  the  above,  a  knowledge  of  £  within  a  interval 
£  +  0.05  was  considered  adequate  for  the  purpose  of  the  present 
investigation . 

Calculations  were  then  carried  out  in  the  neighborhood  of  the 
first  approximation  best  value  of  £  with  the  steps  in  £  which 
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al lowed  us  to  locate  f  with  the  said  approximation. 

The  results  of  the  data  processing  is  summarized  in  Tables  IVa 
and  IVb. 

Tab. IVa  shows  the  best  values  £  and  n  obtained  in  the  case  of 
stiffened  panels  with  countersunk  rivets;  tab. IVb  gives  the  same 
quantities  for  the  stiffened  panels  with  round  head  rivets. 

As  far  as  the  countersunk  rivet  data  is  concerned,  the  scatter 
in  n  is  negligible  and  the  vast  majority  of  values  fall  in  the 
interval  0.105  £  n  $  0.175  and  they  have  been  assigned  the  value 
n  =  0.14. 

For  such  data  it  was  possible  to  obtain  the  cumulative  distri¬ 
bution  of  £.  Log-normal  distribution  of  such  a  variable  is  given 
in  Fig. 9. 

As  far  as  the  data  relevant  to  round  head  rivets  is  concerned, 
a  wider  scatter  in  n  has  been  found  and  the  present  data  does  not 
allow  us  to  drawn  conclusions  in  term  of  the  distribution  of  such 
random  variables. 

Further  data  must  be  collected  to  obtain  a  clear  picture  of  the 
phenomenon. 

3,5  -  LIMITS  OF  THE  PRESENT  ANALYSIS. 

The  results  summarized  in  the  previous  section  have  two  limita¬ 
tions.  The  first  one,  which  is  inherent  in  the  data  processing 
system,  stems  from  the  accuracy  with  which  the  spline  module  works. 

As  already  explained,  the  logic  of  the  present  approach  is  to 
ascribe  any  anomaly  in  the  K-rate  relationship  to  K  evaluation  and 
to  determine  the  best  values  of  flexibility  and  friction  by  minimizing 
such  anomalies.  But,  if  anomalies  in  the  K-rate  relationships  stem 
from  errors  in  da/dn  computation  procedure,  the  best-fit  technique 
tries  to  eliminate  such  anomalies  by  modifying  £  and  n. 

As  a  consequence  of  the  above,  it  is  necessary  to  guarantee  nec[ 
ligible  errors  in  the  output  of  the  spline  module. 

However,  it  is  difficult  to  quantify  the  accuracy  of  the  method 
and  one  must  rely  mainly  on  checks  performed  utilizing  comparison 
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with  known  functions. 

As  explained  in  appendix  1 ,  this  was  the  way  used  in  the 
evaluation  of  the  spline  method. 

Good  results  were  obtained  by  applying  the  spline  module  to 
known  functions  which  present  the  typical  trend  of  a-n  data 
found  in  stiffened  panels,  and  by  comparing  the  spline  derivatives 
with  the  analytical  ones. 

However,  examination  of  the  results  in  terms  of  da/dn  versus 
K  as  in  the  examples  in  Fig. 7,  gave  rise  to  some  doubts  with 
respect  to  the  method  used  to  assess  the  accuracy  of  the  spline 
module . 

It  was  decided  to  proceed  to  a  further  evaluation  based  on  the 
following  methods. 

A  Paris  K-rate  relationship  is  numerically  integrated  utilizing 
SIF  solutions  corresponding  to  stiffened  panels  with  different 
rivet  flexibilities. 

The  a-n  data  so  obtained  is  then  utilized  as  input  to  the  SKESA 
computer  program. 

In  such  a  way,  since  the  actual  value  of  the  flexibility  is 
known  and  since  the  errors  induced  by  the  integration  are  neg¬ 
ligible,  the  difference  between  the  actual  flexibility  and  that 
computed  by  SKESA,  can  be  used  to  quantify  the  error  in  the  spline 
module. 

Fig. 10  shows  the  results  of  such  an  investigation.  As  can  be 
observed,  errors  up  to  10%  in  £  are  to  be  found;  such  errors  are 
always  on  the  same  side  and  the  computed  solution  is  stiffer  than 
the  real  one.  Fig. 11  shows  a  comparison  between  a  typical  da/dn-a 
curve  as  computed  by  the  spline  module  and  the  actual  values  of 
such  a  function. 

The  spline  module  fails  to  reproduce  da/dn  accurately  only  in 
the  two  zones  of  maximum  and  minimum  predicting  a  lower  rate  in 
the  zone  of  the  maximum  and  a  higher  rate  in  the  zone  of  the  min¬ 
imum. 

Improvement  in  the  spline  module  output  might  be  obtained  by 
thickening  the  points  to  be  interpolated  in  the  neighbourhood  of 
the  points  of  extremum. 
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Preliminary  results  of  such  an  approach, of  which  a  typical 
example  can  be  seen  in  Fig. 11,  indicate  that  such  a  way  is 
practicable  even  if  no  well  founded  rule  of  thickening  is  at 
present  available. 

The  second  limitation  stems  from  the  nature  of  the  scatter  in 
the  parameters  £  and  n.  Since  flexibility  and  friction  are  random 
variables,  different  fasteners  in  the  same  stiffened  panel  should 
be  assigned  different  values  of  £  and  n  in  K  computation . Obviously 
such  a  point  of  view  can  hardly  be  implemented  within  the  limits 
of  acceptable  computer  times. 

The  present  approach  which  is  based  on  constant  values  of  £ 
and  n  for  a  given  panel,  gives  on  the  contrary  weighted  averages 
of  such  quantities  which  can  be  effectively  used  to  compute  the 
SIF  and  therefore  the  crack  growth  rate. 

Notwithstanding  such  limitations,  the  approach  implemented  in 
the  SKESA  computer  program  can  be  considered  an  effective  device 
for  analysis  of  crack  growth  data  in  stiffened  panels  in  order  to 
obtain  an  engineering  evaluation  of  the  fastener  flexibility  and 
friction  forces. 

4  -  EVALUATION  OF  EXISTING  METHODS  FOR  PREDICTING  CRACK  GROWTH 
UNDER  VARIABLE  AMPLITUDE  LOADING. 

As  explained  in  section  2 ,  the  next  step  in  the  investigation 
concerns  evaluation  of  existing  methods  for  predicting  crack  growth 
under  variable  amplitude  loading. 

Attention  has  been  centered  mainly  on  the  following  methods, 
namely  the  non-interactive,  Wheeler, | 3 | , Bell-Eidincf f , | 4 | ,  and 
Willenborg, j  2 | methods . 

The  growth  of  a  crack  under  variable  amplitude  loading  is  a 
complicated  phenomenon  not  fully  understood  at  present. 

The  prediction  methods  largely  rely  upon  empiricism  and  their 
ability  to  compare  with  test  data.  In  such  conditions  any  approach 
for  the  evaluation  of  such  methods  must  be  founded  on  comparison 


with  test  data  on  the  basis  of  some  criterium  allowing  us  to 
evaluate  the  fit  between  test  and  computation  data. 

A  first  question  in  this  respect  is  the  fact  that  crack  growth 
is  a  random  phenomenon  whereas  prediction  methods  are  deterministic. 

As  a  consequence,  any  comparison  must  be  judged  on  a  statistical 
basis . 

Further,  all  the  prediction  methods  considered  in  the  present 
investigation  are  cycle  by  cycle  computation  procedures . They  are 
based  on  constant  amplitude  crack  propagation  data  and  account  for 
interaction  effects  by  modifying  the  stress  cycle  at  the  crack  tip 
on  the  basis  of  set  of  rules  which  differ  from  method  to  method  (no 
cycle  modification  for  the  non-interactive  method) . 

Therefore,  if  you  are  interested  in  evaluating  such  rules, you 
need  constant  amplitude  crack  growth  data  which  allows  you  to  make 
an  unbiased  prediction.  In  other  words,  any  error  in  prediction 
needs  to  stem  from  the  rule  and  not  from  bad  constant  amplitude  data. 

As  a  consequence  of  the  above,  the  evaluation  carried  out  in  the 
present  investigation  has  been  based  on: 

-  a  computer  program  which  implements  the  above  mentioned  prediction 
methods.  Such  a  computer  program,  named  CADAV,  has  already  been 
described  in  ref. | 9| ; 

-  a  rationale  to  produce  unbiased  constant  amplitude  crack  growth 
data  and  to  judge  on  a  statistical  basis  the  fit  between  experi¬ 
mental  and  computed  crack  growth  data; 

-  a  set  of  crack  growth  test  data  under  constant  amplitude  loading; 
such  data  has  been  used  to  define  the  constant  of  the  semLempirical 
laws  (Paris, Forman  and  Collipriest)  which  give  analytical  expres 
sion  to  the  K-rate  relationship; 

-  a  set  of  crack  growth  test  data  under  variable  amplitude  loading 
to  compare  with  the  results  of  the  prediction  methods.  Such  data 
has  been  obtained  by  testing  centrally  through  cracked  flat 
specimens  with  standardized  load  sequences  FALSTAFF, | 10 |  and 
TWIST, | 11 | ,  following  the  test  techniques  described  in  appendix 
2. 


L\,l  -  THE  RATIONALE  OF  THE  PRESENT  APPROACH. 


Fig. 12  shows  the  main  steps  of  the  rationale. 

It  comprises  two  main  lines  which  refer  to  constant  amplitude 
data  generation  and  analysis  and  variable  amplitude  data  generation 
and  analysis. 

The  objective  of  the  first  line  is  the  determination  of  the 
semiempirical  laws  (Paris , Forman , Collipriest-Walker)  for  the 
material  used  in  the  variable  amplitude  tests  with  sufficient  ac¬ 
curacy  to  allow  us  to  make  an  unbiased  prediction  of  crack  growth 
under  variable  amplitude  loading. 

To  this  end  starting  from  the  crack  growth  data  obtained  with 
the  constant  amplitude  test  program,  (see  sec. 4. 2),  the  following 
main  operations  have  been  carried  out: 

-  generation  of  K-rate  relationships  in  the  form  given  in  the 
Paris, Forman  and  Collipriest  law  by  means  of  the  SKESA  computer 
program.  SKESA  selects  the  constants  which  define  such  laws 
through  a  regression  analysis  substantiated  by  the  usual  signif 
icance  tests; 

-  computation  of  the  number  of  Nc  cycles  necessary  to  reach  a  given 
crack  length  by  integration  of  the  best  fit  curves  defined  in  the 
previous  step; 

-  comparison  of  the  experimental  crack  growth  data  with  the  data 
computed  by  means  of  a  statistical  analysis  of  the  random  va¬ 
riable  nex/nc  (nex  being  the  test  value  of  nc) .  Such  an  analysis 
consists  in  the  determination  of  the  cumulative  distributions  of 
nex/nc  for  different  crack  lengths  with  data  coming  from  several 
specimens . 

The  first  attempt  at  best-fit  gives  an  unbiased  prediction  only 
if  the  medians  of  the  distributions  corresponding  to  different 
crack  lengths  have,  apart  from  small  differences,  always  the  same 
value  (nex/nc)M  =  1. 

If  this  is  not  so, it  means  that  the  semiempirical  laws  do  not 
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follow  the  same  trend  as  the  test  data  in  some  part  of  the  growth 
interval.  A  typical  cause  of  this  state  of  affairs  is  unprevented 
crack  lip  buckling  in  the  high  crack  length  range  which  increases 
the  crack  growth  rate  with  respect  to  the  unbuckled  condition 
which  the  SIF  solution  refers  to. 

In  these  cases,  some  test  data  must  be  left  out  and  the  steps 
previously  taken  must  be  repeated  on  the  reduced  set  of  a-n  data. 

Once  unbiased  constants  of  the  semiempirical  laws  have  been 
obtained,  evaluation  of  the  variable  amplitude  data  can  be  carried 
out . 

The  main  steps  of  the  variable  amplitude  line  are  the  following: 

-  utilization  of  data  from  damage  growth  tests  in  specimens  subjected, 
to  variable  amplitude  loading. 

These  tests  must  be  conducted  with  specimens  made  of  material 
taken  from  the  same  batch  used  in  constant  amplitude  test  spec¬ 
imens.  Different  standardized  spectra  can  be  used,  namely  FALSTAFF , 
TWIST  and  Gaussian  random;  tests  with  different  values  of  refer¬ 
ence  stress  of  spectrum  can  be  performed.  For  each  case  -assigned 
specimen  and  spectrum  at  a  given  value  of  reference  stress-  a 
sufficient  number  of  tests  are  performed  to  provide  a  meaningful 
statistical  analysis; 

-  prediction  of  damage  growth  with  the  non-interactive,  Wheeler, 
Bell-Eidinof f  and  Willenborg  methods  for  each  test  case. 

The  computer  program  CADAV  is  used  for  this  purpose; 

-  comparison  of  the  experimental  damage  growth  data  with  the  theo¬ 
retical  data.  The  comparison  is  obtained  by  determining  the  cumu 
lative  distribution  of  the  random  variable  Fex/Fc,  Fex  being  the 
number  of  flights  in  which  the  damage  grows  to  an  assigned  value 
of  size  and  Fc  being  the  same  data  obtained  with  the  prediction 
method. Such  a  distribution  is  obtained  for  different  intervals  of 
damage  propagation,  utilizing  an  adequate  number  of  tests.  The 
median  value  of  Fex/Fc  can  be  used  to  measure  the  fit  between 
test  data  and  prediction  methods.  Median  values  near  one  for  all 
the  distributions  indicate  an  excellent  capability  of  the  method 
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to  predict  crack  growth  under  variable  amplitude  loading. 

Scatter  around  the  median  value  is  another  important  piece 
of  information  that  can  be  obtained  from  inspection  of  the  dis 
tributions.  In  particular,  comparison  of  scatter  between  constant 
and  variable  amplitude  data  is  of  fundamental  importance. 

4,  2  -  EXPERIMENTAL  RESULTS  AND  ANALYSIS. 

4,2.1  -  Constant  amplitude. 


Following  the  rationale  of  the  previous  section,  constant  ampli 
tude  tests  have  been  performed  in  order  to  obtain  the  expression  of  sc 
iempirical  laws. 

Twelve  flat  specimens  with  central  through  crack  made  in  2024-T3 
have  been  tested  for  various  values  of  the  ratio  R=Smin/SmaX*  The 
test  apparatus  and  test  procedures  are  shown  in  appendix  2. 

A  detailed  explanation  of  data  processing  is  given  in  ref.|9|. 

Here  the  final  results  are  reported  as  far  as  the  Forman  law 
determination  is  concerned.  Fig. 13  shows  data  for  all  the  tests  and 
the  best-fit  straight  line  as  obtained  by  regression  analysis, 
together  with  the  scatter  band  of  10%-90%.  In  the  same  figure,  the 
values  of  the  typical  constants  of  the  law  are  also  shown. 

Fig. 14  shows  the  best-fit  straight  line  relevant  for  tests  with 
different  values  of  parameter  R.  The  three  straight  lines  lie  very 
close  to  each  other  and  demonstrate  the  effectiveness  of  Forman's 
law  to  describe  damage  growth  phenomenon  at  different  values  of  R; 
this  is  an  important  requirement  for  a  law  which  must  be  used  in 
prediction  methods  for  variable  amplitude  loading. 

The  Forman's  law  so  obtained  has  been  further  analyzed  with 
respect  to  its  capability  to  predict  the  number  of  cycles  needed 
to  propagate  the  damage  from  the  initial  dimension  to  an  assigned 
dimension.  The  results  are  given  in  Fig. 15a  which  shows  the  dis¬ 
tribution  on  normal  probability  paper  of  the  variable  log(nex/nc) 
for  two  intervals  of  cracks  growth. 
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For  the  propagation  from  aQ  =  6  mm.  to  a  =  30  mm. ,  the  distri 
bution  is  normal  with  significant  accuracy,  and  the  median  value 
is  very  close  to  zero;  it  means  that  the  Forman's  law  obtained 
is  fairly  representative  of  the  crack  growth  phenomenon  in  this 
crack  interval. 

For  the  propagation  from  aQ  =  6  mm.  to  a  1 8  mm.,  the  law  does 
not  work  so  accurately  as  in  the  previous  case. 

This  state  of  affairs  may  be  ascribed  to  higher  relative  error 
which  affects  the  test  data  in  the  low  crack  length  range  since 
the  error  in  measuring  crack  dimension  is  constant.  Nevertheless, 
the  prediction  can  still  be  considered  satisfactory  and  the  law 
so  determined  has  been  used  in  the  prediction  methods  for  crack 
propagation  under  variable  amplitude  loading. 

4.2.2  ~  Variable  amplitude  loading. 

For  the  second  aspect  of  this  research,  seven  tests  have  been 
performed  using  FALSTAFF  spectrum, | 10 | , ^  with  Sj^x  =  2^5  MPa  on 
flat  panels  centrally  through  cracked  made  with  material  2024-T3, 
drawn  from  the  same  batch  used  for  specimens  of  constant  amplitude 
tests.  The  test  apparatus  and  procedure  are  described  in  appendix 
2. 

The  results  of  the  analysis  are  shown  in  Fig. 16  and  Fig. 17. 

Fig. 16  shows  the  crack  length  versus  the  number  of  flights  as 
obtained  in  the  tests  of  the  seven  specimens  and  as  predicted  by 
the  non-interactive,  Wheeler,  Bell-Eidinof f  and  Willenborg  meth¬ 
ods. 

Fig. 17  shows  the  distribution  on  normal  probability  paper  of 
the  variable  log(Fex/Fc)  for  two  intervals  of  crack  damage. 

Such  a  variable  conforms  to  a  normal  distribution  for  the  two 
intervals.  The  scatter  of  phenomena,  measured  from  the  standard 
deviation  of  the  variable,  is  pratically  the  same  as  that  found 

(1)  Preliminary  testa  using  the  TWIST  spectrum  have  been  also  per 
formed.  The  results  are  given  in  ref.\9\. 
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in  constant  amplitude  data.  Such  a  conclusion  seems  to  indicate 
no  particular  effects  of  delay  or  acceleration  effects  due  to 
peak  load  on  the  scatter. 

For  each  prediction  method  the  median  of  the  distribution  (50% 
value)  provides  a  good  indication  of  the  correctness  of  the  meth¬ 
od. 

The  results  are  still  limited  and  therefore  no  general  conclu 
sions  can  be  drawn,  but  certain  trends  in  the  flat  panels  are 
well  established.  The  non-interactive  method  predicts  life  in  the 
safe  side  with  a  factor  of  about  6  in  accordance  with  the  pronounced 
delay  effects  due  to  the  high  peak  loads  typical  of  the  FALSTAFF 
spectrum. 

The  Wheeler  method,  with  an  adequate  selection  of  the  value  of 
its  plastic  zone  characteristics  constant,  m  =  1.9,  produces  an 
"exact"  prediction,  namely  the  mean  value  of  random  variable  log 
(Fex/Fc)  equal  to  zero. 

The  Willenborg  method  gives  unconservative  results  with  the 
FALSTAFF  spectrum,  where  the  presence  of  high  compressive  loads 
produces  crack  acceleration  which  is  not  taken  into  account  by 
this  method. 

The  Bell-Eidinof f  method  does  not  tally  satisfactorily  with  ex 
perimental  data.  A  possible  cause  of  this  state  of  affairs  may  be 
ascribed  to  the  fact  that  the  values  of  the  empirical  constants 
which  are  taken  from  the  existing  literature  and  which  are  used 
to  characterize  the  delay  and  acceleration  effects,  might  not  work 
very  well  with  the  material  used  in  the  present  investigation.  The 
values  of  such  a  constant  directly  determined  by  ad-hoc  tests  might 
improve  the  prediction. 

The  test  program  on  stiffened  panels  under  standardized  spectra 
is  now  in  progress  and  in  this  case  the  experimental  data  obtained 
so  far  shows  a  larger  scatter,  Fig.  18,  with  respect  to  the  flat 
panels. 

The  data  of  Fig.  18  is  still  preliminary  and  no  comparison  with 
the  calculation  method  has  yet  been  made. 
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5  -  CONCLUSIONS 


An  investigation  on  fatigue  crack  growth  in  aircraft  structures 
is  currently  being  carried  out  at  the  Institute  of  Aeronautics  of 
the  University  of  Pisa,  Italy,  under  contract  DA-ERO-78-G-1 07 . 

The  final  goal  of  the  investigation  is  the  development  of  a 
reliable  method  to  predict  the  growth  of  a  crack  in  built-up 
structures  under  operational  load  conditions. 

The  first  year's  results,  presented  in  this  paper,  mainly 
concern  the  two  following  topics,  namely 

-  crack  growth  in  riveted  stiffened  panels  under  constant  amplitude 
loading 

-  crack  growth  in  plain  sheets  under  variable  amplitude  loading. 

As  far  as  the  first  topic  is  concerned,  the  results  are  relevant 
mainly  to  the  evaluation  of  the  rivet  flexibility  and  friction 
forces  between  stringer  and  sheet-cover.  The  results  of  crack 
propagation  tests  of  stiffened  panels  part  of  which  are  riveted 
with  countersunk  rivets,  the  others  with  round  heat  rivets,  have 
been  analysed  with  a  SKESA  computer  program  developed  ad-hoc, which 
allows  us  to  obtain  the  average  values  of  the  rivet  flexibility 
and  friction  forces.  SKESA  computes  such  quantities  by  searching 
for  the  best  regression  in  correlating  the  da/dn  crack  rate  with 
AK. 

It  works  with  a  good  degree  of  accuracy, the  errors  in  final 
output  being  generally  of  the  order  of  5%  and  in  each  case  not  higher 
than  10%.  They  depend  mainly  on  the  spline  module  which  computes 
da/dn  from  the  test  data  a-n,  and  could  be  reduced  by  a  rational 
thickening  of  the  points  (a,n)  to  be  interpolated  in  the  neigh¬ 
bourhood  of  the  flex  points  in  the  a-n  curves.  This  state  of  af¬ 
fairs  implies  a  modification  in  the  test  technique  with  a  different 
selection  of  the  reading  spacings  in  the  measurement  of  the  crack 
length. 

Notwithstanding  such  errors  the  SKESA  computer  program  represents 
a  valuable  device  which  allows  us  to  obtain  quantitative  information 
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on  flexibility  and  friction  which  can  be  then  confidently  used  for 
analysis  purpose. 

In  particular,  it  has  already  been  possible  to  obtain  conclusive 
results  as  far  as  countersunk  rivets  are  concerned,  by  defining 
the  cumulative  distribution  of  the  flexibility  parameter  £ ,  and 
finding  the  range  in  which  the  friction  parameter  n  lies.  On  the 
contrary,  in  the  case  of  round  head  rivets,  the  pattern  is  more 
scattered  and  there  is  still  not  enough  data  available  at  present 
for  statistical  analysis. 

New  research  in  this  area  will  be  carried  out,  in  the  continua 
tion  of  the  program,  by  improving,  if  possible,  the  spline  module 
and  by  producing  further  test  data  with  stiffened  panels  construc¬ 
ted  with  different  types  of  fasteners. 

As  far  as  the  second  topic  is  concerned,  four  methods  for  com¬ 
puting  crack  growth  under  variable  amplitude  loading  have  been 
evaluated. 

Such  methods  are  the  non-interactive,  Wheeler,  Bell-Eidinof f  and 
Willenborg  methods. 

Evaluation  has  been  carried  out  firstly  by  producing  constant 
amplitude  crack  growth  data  to  obtain  unbiased  semiempirical  laws 
(of  the  Paris,  Forman  and  Collipriest  type)  and  then  performing 
variable  amplitude  crack  growth  tests  with  the  FALSTAFF  spectrum 
on  specimens  made  with  material  coming  from  the  same  batch  as  the 
constant  amplitude  test  specimens.  The  numbers  of  flights  necessary 
to  increase  the  crack  in  a  given  interval  Fex  as  obtained  by 
testing  different  specimens,  are  compared  with  predictions  Fc  as 
obtained  from  the  CADAV  computer  program,  which  implements  the 
above  mentioned  method  by  utilizing  the  previously  determined 
semiempirical  laws.  The  comparison  is  made  on  a  statistical  basis 
by  finding  the  distribution  of  Log  Fex/Fc  for  the  different  meth¬ 
ods  and  for  different  crack  growth  ranges.  The  median  of  the 
distribution  is  then  used  as  measurement  of  the  accuracy  of  each 
method . 

The  results  obtained  up  to  now  indicate  that  the  Wheeler  meth¬ 
od  can  be  confidently  used  at  least  with  the  FALSTAFF  spectrum. 
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An  important  conclusion  of  this  analysis  concerns  the  scatter 
which  has  been  found  to  be  of  the  same  order  both  in  constant 
and  variable  amplitude  tests. 

Further  research  in  this  area  will  be  carried  out  by  per¬ 
forming  tests  both  on  plain  specimens  and  stiffened  panels  util¬ 
izing  new  loading  spectra  (TWIST  and  Gaussian  random  loads)  and 
increasing  the  test  data  obtained  with  the  FALSTAFF  spectrum. 
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APPENDIX  1  -  SPLINE  TECHNIQUE 


One  of  the  most  difficult  problems  in  analysing  crack  propa¬ 
gation  data  is  the  definition  of  the  crack  rate,  that  is  how  to 
"differentiate"  the  discrete  raw  data  obtained  in  experiments. 

The  problem  is  of  great  importance,  since  the  parameters  of  the 
semi-empirical  laws  of  crack  propagation,  which  link  the  crack 
rate  da/dn  with  the  AK  stress  intensity  factor  range,  are 
strongly  influenced  by  the  accuracy  with  which  the  derivatives 
da/dn  are  computed. 

In  the  past  a  common  solution  for  the  problem,  especially  for 
data  obtained  from  simple  sheet  panels,  lying  in  this  case  on  a 
smooth  curve,  has  been  provided  by  techniques  such  as  finite 
difference  methods.  But  the  data  points  are  subjected  to  a  degree 
of  scatter  which  results  in  large  errors  of  growth  rate  when 
using  this  method.  So  the  technique  of  fitting  a  continuous  curve 
to  the  data  points  using  the  least  squares  concept  has  been 
developed,  and  the  crack  growth  rate  can  be  obtained  simply  by 
differentiating  the  analytical  expression  of  a  (n) . 

At  first,  a  method  of  fitting  a  cubic  polynomial  function  was 
adopted  but,  for  the  case  of  stiffened  panels, it  was  decided  to 
make  use  of  spline  function  (piece-wise  continuous  polynomial) 
and  a  very  good  fitting  was  obtained. 

The  method  adopted,  explained  in  detail  in  |12|,  consists  in 
fitting  in  the  least  square  sense  a  curve  of  the  type: 


n  =  X  -  Y  ln(a) 


Z  n 

|  +  *  *p(a-Kp) 
P=i 


(1) 


to  the  experimental  data  points  (aj,nj).  The  first  part  of  the 
expression  is  a  function  of  a  class  suggested  by  a  cumulative 
concept  of  damage,  while  the  second  is  the  spline  function, with 
polynomials  of  degree  l  connected  at  knots  Kp.  Following  the 
usual  method  of  minimizing  the  sum  of  squares  of  residuals,  a 
computer  program  carries  out  the  computations  for  selecting  the 
best-fit  values  of  the  parameters  X,Y,Z,Xp  of  eq.(1). 
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Tests  have  been  performed  to  choose  proper  values  for  the 
number  of  knots  and  for  the  degree  of  the  polynomials  and  the 
results  show  that  the  values  selected  of  4  knots  for  long  cracks 
and  3  for  short  ones,  together  with  4  for  the  exponent,  lead  to 
satisfactorily  small  errors. 

Tests  have  also  been  carried  out  by  monitoring  the  errors 
made  when  comparing  the  derivatives  of  known  functions  (such 
as  polynomials,  trigonometric  functions,  etc..)  with  their 
approximations  obtained  by  using  the  spline  function  method. 

The  results  have  been  considered  satisfactory  and  so  the  meth¬ 
od  has  been  adopted  as  a  standard  for  getting  crack  growth  rates 
from  experimental  points. 

Other  comparisons  have  also  been  made  with  functions  represent 
ing  propagation  curves  of  ideal  stiffened  panels  under  constant 
amplitude  fatigue  loading;  the  results  have  been  discussed  in  the 
text,  paragraph  3.5. 
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APPENDIX  2  -  TEST  APPARATUS 

The  equipment  prepared  for  the  damage  propagation  tests  is  com 
posed  of  two  different  load  machines,  of  the  load  spectrum  gener¬ 
ation  system  (in  variable  amplitude  loading  tests)  and  the  defect 
growth  measuring  system. 

The  first  piece  of  loading  equipment  is  composed  of  a  Servotest 
hydraulic  actuator,  controlled  in  closed  loop  through  a  set  of  six 
servovalves,  capable  of  +  500  KN  or  +  100  KN  (static  and  dynamic) , 
fed  by  a  125  HP  power  pack  with  an  oil  flow  rate  of  190  Lt./min. 

For  constant  amplitude  tests,  the  control  unit  of  the  machine 
directly  drives  the  execution  of  load;  for  variable  amplitude 
loading  spectra  an  external  input  is  available.  The  actuator  is 
mounted  on  a  very  stiff  rig  and  the  specimen  is  clamped  to  the  rig 
and  to  the  load  cell  of  the  actuator  with  grips  which  hold  the 
specimen  by  friction.  Moreover,  anti-buckling  guides  are  available 
both  for  flat  specimens  and  stiffened  specimens. 

The  second  piece  of  equipment  is  composed  of  a  Servotest  hydrau 
lie  actuator,  controlled  in  closed  loop  through  a  set  of  two 
servovalves,  capable  of  +  250  KN  or  +  50  KN  (static  and  dynamic) . 
This  actuator  is  in  every  way  similar  to  the  first,  is  fed  from 
the  same  power  pack  and  is  mounted  on  a  similar  rig. 

The  electric  signal,  which  drives  the  actuators  in  the  case  of 
variable  amplitude  loading,  is  generated  from  a  process  computer 
PDP  11/34.  Preliminary  tests  for  determining  the  maximum  working 
frequency  are  performed  for  each  condition,  loading  spectrum  and 
flexibility  of  specimen. 

The  input  signal  in  the  actuator,  representing  the  desired 
spectrum,  is  compared  by  PDP  11/34  with  the  response  of  the  ac¬ 
tuator,  namely  the  output  of  the  load  cell.  The  working  frequency 
selected  is  the  highest  one  for  which  the  difference  between 
input  and  output  is  not  significant.  During  the  tests  the  output 
of  load  cell  is  also  continuously  fed  to  the  PDP  11/34.  In  this 
way,  checks  of  coincidence  of  the  two  signals  can  be  periodically 
performed.  Moreover,  if  the  machine  stops  because  the  maximum  and 
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minimum  load  limits  are  surpassed  or  for  other  reasons,  the  PDP 
11/34  commands  the  arrest  of  the  input  signal  and  prints  the 
number  of  flights  at  that  moment. 

The  defect  growth  measuring  system  is  formed  essentially  of 
a  30  magnification  telescope,  whose  position  can  be  regulated, 
with  a  movement  parallel  to  the  propagation  direction,  along 
a  bar,  which  has  a  millimetric  scale;  a  vernier  is  applied  to 
the  telescope  bearing  block,  allowing  the  defect  dimension  measure 
to  be  taken  with  a  0.1  mm.  precision.  The  number  of  cycles,  or 
more  exactly  the  number  of  flights  in  the  "flight-by-flight" 
spectra  case,  at  which  the  defect  has  reached  the  measured 
dimension,  is  written  on  a  printer  connected  to  a  counter  in 
which  the  actuator  load  cell  signal  is  fed:  the  printer  is  ac¬ 
tivated  directly  by  the  operator  who  makes  the  observations.  In 
the  case  in  which  the  spectrum  is  generated  by  the  process  com¬ 
puter,  the  flight  number  is  directly  printed  on  the  telescope  type 
when  the  operator  sends  a  proper  signal  to  the  process  computer 
input.  The  whole  measuring  system  allows  us  to  take  highly 
accurate  measurements,  certainly  adequate  for  the  purposes  of 
this  research. 
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Fig.1  -  Typical  anomalies  in  K-rate  relationships  in  stiffened  panels 
due  to  errors  both  in  SXF  and  da/dn  computation. 
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fASItNtK  mXIBILIIY  —  IHtORUICAL  APPROACH 
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Fig. 5b  -  Crack  growth  data  for  stiffened  panels  riveted  with  round  head  rivets.  New  series  tests. 
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Fig. 7c  -  Test  CR1C  -  1:  da/dn  vs.  AK  for  different  values  of 

£  and  n,  plotted  as  explained  in  Fig. 7a.  In  this  case 
loop-type  anomalies  are  present  also  at  the  higher 
values  of  FDATA. 
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g.9  -  Log  C  distribution.  The  curve  A  is  rele 

vant  to  n  =  0.14.  The  median  value  is  £:= 

=  0.49  (Ci  =  0.686)  and  the  standard  de 
viation  is  02=0.219  (oj  =0.307).  The  re¬ 
sults  indicate  that  the  flexibilities  found 
with  this  approach  are  noteworthy  lower  than 
those  predicted  by  Fig.  3  or  Douglas  formula. 
The  curve  B  has  been  obtained  considering  all 
the  l  belonging  to  the  same  population  irre¬ 
spective  of  the  values  of  n  . 


Fig. 10  -  Errors  found  in  C  evaluation  due  to  errors  in  the 
spline  module.  Within  range  of  pratical  interest 
0.2se2S°*8  the  absolute  error  is  lower  than  0.06 
and  the  relative  error  is  lower  than  12%. 
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Fig. 11  -  Typical  errors  of  the  spline  interpolation 

technique.  The  curve  with  equidistant  points 
is  representative  of  the  present  approach. 

The  curve  with  thickened  points  gives  an  idea 
of  potential  improvements  that  can  be  obtained 
by  a  modified  test  technique  which  provides  more 
information  in  the  neighbourhood  of  maximum  and 
minimum  points. 
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Fig. 13  -  Regression  analysis  on  experimental 
data  for  evaluating  Forman's  law. 
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Fig. 14  -  Regression  analysis  on  experimental 
data  for  evaluating  Forman's  law. 
Best-fit  straigth  lines  for  data 
relevant  to  R=cost. 
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Tab. Ill  -  Results  of  new  series  tests. 
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52.8 

53525. 

54.5 

55475. 

57.0 

58390. 

58.6 

60100. 

60.8 

62460. 

61.7 

63610. 

63.2 

65165. 

65.6 

68065 . 

66.8 

69605. 

68.8 

72145. 

70.4 

74125. 

73.1 

7(7425. 

75.0 

79685. 

78.9 

6^4100 . 

81.5 

86840. 

83.4 

£8835. 

85.4 

90770. 

87.8 

93140. 

90.6 

95795. 

m&Wmk 

97415. 

100425. 

98.4 

102660. 

100.8 

104675. 

103.6 

106965. 

106.4 

109075. 

109.4 

111280. 

111.5 

112750. 

CRIC-5 

a 

n 

16.1 

7930. 

16.5 

11950. 

17.5 

19985. 

18.1 

25390. 

19.1 

31600. 

19.7 

34920. 

20 . 1 

37885. 

21.0 

22.1 

49070. 

23.0 

55000. 

23.8 

58460 . 

25 . 3 

66020 . 

25.9 

69415. 

27.4 

76240. 

28.1 

79135. 

29.5 

85440. 

30.5 

88990. 

32.2 

96010. 

33.6 

101605. 

35.1 

106700. 

37.5 

115010. 

38.6 

1 18535. 

41.7 

127960. 

44.3 

135295. 

46.5 

141605. 

49.1 

148595. 

49.8 

150895. 

52.9 

159065. 

167390. 

El  ts 

173900. 

60.5 

180450. 

62.8 

186890. 

64.5 

192625. 

65.8 

196335. 

68.8 

205510. 

70.5 

210720. 

73.1 

218535. 

75.0 

223920. 

77.9 

232495. 

79.3 

236495. 

82.6 

245590. 

84.3 

249965. 

86.2 

255000. 

88.9 

261935, 

91.2 

267085. 

cure 

-6 

a 

n 

16.4 

17570. 

17.1 

21305. 

18.1 

29670. 

19.5 

38845. 

20.7 

49420. 

21.5 

53800. 

22.8 

60360. 

23.5 

660 1 0 . 

24.3 

70120. 

25.7 

77070. 

27.5 

85500 . 

29.7 

94450. 

30.5 

98555 . 

31.5 

102675. 

32.5 

106825. 

33.6 

111260. 

34.5 

115705. 

36  ♦  0 

120230. 

36. 8 

123405. 

38.2 

128000. 

39.3 

131640. 

40.6 

136750. 

43.0 

143500. 

44.9 

149650. 

46.2 

153070, 

47.3 

156160. 

49.4 

161650. 

50.4 

164890. 

52.2 

169065. 

54.5 

175105. 

55.5 

177775. 

180840. 

58.4 

184490. 

59.5 

187805. 

60.5 

190655. 

61.2 

192200. 

196485. 

64.1 

199945. 

65.3 

203185. 

66.5 

208520. 

68 » 6 

212015. 

70.4 

216985. 

71.3 

219000. 

73.0 

223625. 

74.5 

227380. 

« 
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CHIC-9 

a 

n 

16.4 

7430. 

17.6 

10850. 

IB. 8 

13035. 

19.3 

15820. 

20.4 

18235. 

22.0 

22005. 

23.0 

23675. 

24.1 

25950. 

25.9 

28675. 

27.4 

31760. 

28.4 

33030. 

29.3 

35260. 

30.4 

37305. 

31.2 

38910. 

32.7 

41580. 

34 . 0 

43635. 

35.3 

45895. 

36.8 

48215. 

39.3 

51610. 

41.0 

54325. 

42.8 

56780. 

44.4 

58880. 

46.0 

60995. 

47.5 

62835. 

51.7 

67740. 

53.5 

69825. 

55.0 

71785. 

57.2 

74490. 

58.8 

76045. 

60.3 

77980. 

62.2 

80130. 

63.7 

82145. 

65. 1 

84125. 

66  4  3 

85590. 

68.1 

88065. 

69.7 

90265. 

70.8 

91775. 

71.9 

93290. 

73.4 

95380. 

75.3 

97675. 

76.8 

99585. 

78.4 

101575. 

81.1 

104865. 

82.5 

106450. 

84.7 

108980. 

cure- 7 

a 

n 

1  8 . 3 

7100. 

1  9 . 5 

11000. 

20.7 

13730. 

22.0 

16870. 

24 . 0 

20845. 

26.0 

25820 . 

28.0 

29515. 

29.4 

32830. 

32 . 0 

36400. 

33 . 5 

38885. 

34 . 5 

40695. 

35 . 8 

42540. 

37.0 

44705. 

39.1 

47030. 

40.6 

49090. 

42.0 

50885. 

44.0 

53350 . 

46.7 

56695. 

48.7 

59160. 

50.0 

60780. 

52.0 

63020. 

53.7 

64640. 

54.9 

66005. 

56.2 

67400. 

57.2 

68600. 

58.1 

69810. 

59.2 

71255. 

60.5 

72495. 

62.0 

74645. 

63.3 

77055. 

64.8 

78255. 

65.7 

79660. 

66.8 

80935. 

67.8 

82160. 

69.1 

83460. 

70.5 

86195. 

72.4 

88510. 

75.6 

92455. 

78.2 

95625. 

80.7 

98570. 

83.0 

101125. 

85.4 

103710. 

87.3 

105720. 

88.9 

107280. 

90,6 

108890. 

cruc-a 

a 

n 

23 . 3 

6765. 

24 . 5 

9210. 

25 . 6 

11430. 

26 . 5 

1 3215. 

27.9 

16115. 

28.5 

17245. 

29.8 

19730. 

32 . 5 

24915. 

33.2 

26135. 

34.8 

28105. 

36.3 

30700, 

37.5 

32 1 80 . 

40 . 5 

36435. 

43.5 

40615. 

44.9 

42320. 

46.0 

43705. 

49.2 

47620. 

50 . 5 

49300. 

52.6 

52850 . 

56.4 

56275. 

59.0 

59255 . 

60.3 

60235. 

62.4 

62845. 

63.2 

64160. 

64.3 

65705. 

66  ♦  3 

68260. 

67,8 

70205. 

70.5 

73665. 

71,7 

75160, 

74.3 

78335. 

75.2 

79430. 

77.0 

81715. 

78.7 

83595. 

81.8 

86790. 

83.2 

88890. 

84.3 

91025. 

86.8 

92180. 

91.8 

95735, 

93.2 

98755. 

94.3 

99855. 

97.1 

101710. 

102.6 

105820. 

106.1 

108570. 

CKIC-IO 


CHIC- 1 1 


CHIC-1  2 


a 

n 

17.6 

13410. 

18.5 

18820. 

19.4 

26740. 

20.6 

34100. 

22.0 

41500. 

22.8 

47320. 

24.1 

53610, 

25.3 

59420. 

26.8 

66390. 

27.8 

73780. 

30.3 

82500. 

31.7 

89310. 

34.0 

96650. 

36.8 

103090. 

37.5 

108510. 

39.8 

115430. 

42.0 

121020. 

43.1 

125580, 

46.2 

139070. 

48.8 

140180. 

51.8 

148810, 

54.0 

155200, 

56.7 

162450. 

58.4 

1 66630 ♦ 

60.3 

171380. 

62.1 

176210. 

63.6 

180770. 

ran 

185040. 

67.3 

190070. 

69.8 

197030. 

71.0 

200560. 

205540. 

74.8 

210390. 

76.2 

214000. 

77.7 

217820. 

79.6 

221610. 

81.4 

226180. 

83.1 

229740. 

84.9 

233520. 

86.9 

238500. 

89.8 

244490. 

250110. 

94.7 

254880. 

97.2 

259510. 

99.8 

263590. 

28.0 
28 .  6 
9. 
JO. 3 

31 . 4 

32 . 8 

34.4 

36.5 
38.3 

39.8 
41 . 1 


43.9 

46.2 
47.5 

49.3 


53.8 

56.0 


61.3 

63 . 3 
3 

66.9 
69.6 
71.2 

73.9 


76.6 

77.9 

78.8 

79.9 
82.0 

83.3 
85.0 

86.3 

87.3 


7645. 
1  I  700 . 
17115. 
18745. 
21855. 
24670. 
26705. 
28135. 
31285. 
33000. 
35100. 
38160. 
41095. 
44705. 
47610. 
50140. 
52105. 
5 

56480. 

59565. 

61595. 

63890. 

67640. 

70010. 

72585. 

74855. 

77475. 

79880. 

82725, 

85505. 

87575. 

91420. 

94750. 

96975. 

99505. 

100895. 

102715. 

103865. 

105195. 

107880. 

109680. 

111230. 

112830. 

113920. 


18.8 

20.5 

21.7 
23.0 

24.8 

26.9 

28.3 

29.6 

30.9 

32.3 

33.7 

35.4 

36.5 

37.8 
39.0 
41.2 

43.1 

44.2 

45.7 

47.7 

49.7 

51.7 


58.7 
60.2 

61.8 

63.3 
64.8 
65 

67.3 

68.3 

69.3 
71.0 

71.7 

73.4 

74.1 

75.8 

77.2 

78.2 
79.1 

80.8 
82.8 
83.8 


n 

7000. 

10700. 

13485. 

16135. 

19670. 

22910. 

25190. 

27365. 

29010. 

31175. 

33380. 

35860. 

37390. 

39200. 

40805. 

43715. 

45980. 

47385. 

49160. 

51570. 

53590. 

55890. 

57855. 

59945. 


63665. 

65505. 

67000. 

69165. 

71080. 

72590. 

74695. 

76260. 

78045. 

80775. 

83245. 

84615. 

85810. 

87505. 

89525. 

90920. 

92250. 

94550. 

96220. 

97265. 
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Okie- 

-13 

cuic- 

-14 

CRIC-1 5 

a 

n 

a 

n 

a 

n 

IB.  1 

10500. 

1 6 . 5 

6915. 

16.5 

10715. 

19.9 

15730. 

17.5 

i 0695. 

16.9 

13875. 

1*1?  •  2 

21995. 

18.9 

14650. 

17.5 

17140. 

24.2 

25915. 

20.4 

18830. 

17.9 

19455. 

25 . 8 

29075. 

21.6 

21945. 

18.5 

21890. 

27. 1 

31285. 

22.5 

24140. 

19.5 

25745. 

28.4 

33420. 

25.6 

31375. 

20.2 

29570. 

30 . 4 

36945. 

27.1 

34595. 

20.9 

32660. 

31.8 

38990. 

28.2 

36605. 

21.9 

35565. 

34.5 

42875. 

29.5 

39145, 

22.7 

38095. 

36.2 

45015. 

30.7 

41230. 

24.0 

37.6 

46525. 

31 . 7 

43150. 

25.3 

46910. 

39.8 

49015. 

33.0 

45265. 

26.5 

50915. 

42.5 

51855. 

33.8 

46560. 

27.5 

54995. 

43.5 

53090. 

34.5 

47860. 

28.7 

58595. 

55175. 

36.8 

i 

61625. 

47.7 

57355. 

37.7 

52715. 

30.7 

64615. 

49.6 

59435. 

39.3 

55005. 

32.3 

67650. 

50.9 

60805. 

40.7 

56890. 

34.0 

70970. 

62635. 

42.5 

59300. 

35.6 

73715. 

64775. 

45.3 

63000. 

37.8 

76895. 

E" 

66395. 

48.0 

66490. 

41.3 

81130. 

58.5 

67690. 

49.7 

68500. 

43.3 

83965. 

60.4 

69570. 

51 . 1 

70140. 

45.8 

86555. 

61.4 

70570. 

52.9 

72280. 

49.6 

90810. 

A3  ♦  3 

72275. 

54.3 

73820. 

93665. 

65.0 

73805. 

55.9 

75720. 

54.8 

95895. 

66.6 

75350. 

58.2 

78475. 

ITTV 

97715. 

68.8 

77665. 

59.5 

80190. 

57.8 

98875. 

69.7 

78490. 

61.3 

82500. 

100240. 

70.7 

79415. 

63.6 

85425. 

60.4 

101230. 

71.9 

80520. 

65.2 

87200. 

61.6 

102370. 

72.8 

81505. 

66.1 

88720. 

62.7 

103365. 

74.2 

82905. 

67.6 

90450. 

64.1 

104695, 

84065. 

70.1 

94130. 

65.1 

105600, 

76.8 

85220. 

71.5 

95880. 

66.9 

107550. 

78.8 

86675. 

73.7 

98705. 

70.0 

110770. 

80.4 

88750. 

100940. 

72.3 

113260. 

83.7 

91765. 

76.1 

101870, 

74.3 

117215. 

85.9 

93715. 

78.4 

119800. 

89.4 

96920. 

79,7 

120945. 

97.9 

103755. 

80.5 

121685. 

■  1 

81,2 

122285 . 

83.2 

124045. 

Tab. Ill 
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5 


TEST 

I2 

CRIC-4 

0.7 

CRIC-8 

0.6 

CRIC-10 

0.75 

CRIC-2 

0.45 

CRIC-3 

0.5 

GSF-13 

0.6 

ALP-6 

0.  5 

GFP-2 

0.  9 

ALP- 4 

0.5 

ALP-7 

0.9 

GSF-5 

0.45 

TEST 

l2 

CRIC-9 

0.45 

CRIC-5 

0.3 

CRIC-6 

1  .0 

CRIC-I 

0.3 

CRIC-7 

0.4 

GFP-3 

0.3 

ALP-5 

0.25 

ALP -3 

0.65 

GSF-1 2 

0.3 

GSF-2 

0.4 

GFP-1 

0.95  , 

0.14 


0.14 


0.14 


0.14 


0.14 


0.14 


0.21 


0.14 


0.14 


0.07 


0.07 


£2  =  Flex. /Flex.  Douglas 
n  -  Pf/Ar>Sy  Sy  =  36  Kg/nun2 


Tab.IVa  -  Optimum  values  of  e;  and  n  parameters  for  countersunk 
rivets  tests. 


TEST 

1  ' 

h 

n 

TEST 

h 

n 

GSF-1 

0.65 

0.1  4 

GSF-4 

0.08 

0.14 

GSF-3 

0.2 

0.14 

GSF-7 

0.2 

0.21 

ALP-2 

0.6 

0.07 

GFP-8 

0.05 

0.07 

ALP-1 

0.2 

0.21 

GSF-1 1 

0.5 

0.21 

GSF-8 

0.7 

0.07 

GSF-9 

0.02 

0.07 

GSF-1 0 

0.22 

0.21 

CRIC-1 1 

0.6 

0.1  4 

CRIC-1 2 

0. 

0.21 

CRIC-1 3 

0.6 

0.07 

CRIC-14 

0.4 

0.1  4 

CRIC-1 5 

0.5 

0.07 

CRIC-1 6 

0.7 

0.07 

=  Flex. /Flex.  Douglas 
n  =  Pf/Ar*Sy  Sy  =  36  Kg/mm2 


Tab.IVb  -  Optimum  values  of  £  and  n  parameters  for  round  head 
rivets  tests. 
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